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ABSTRACT 

Recent direct imaging discoveries suggest a new class of massive, distant planets around A stars. 
These widely separated giants have been interpreted as signs of planet formation driven by gravita- 
tional instability, but the viability of this mechanism is not clear cut. In this paper, we first discuss the 
local requirements for fragmentation and the initial fragment mass scales. We then consider whether 
the fragment's subsequent growth can be terminated within the planetary mass regime. Finally, we 
place disks in the larger context of star formation and disk evolution models. We find that in order for 
gravitational instability to produce planets, disks must be atypically cold in order to reduce the initial 
fragment mass. In addition, fragmentation must occur during a narrow window of disk evolution, 
after infall has mostly ceased, but while the disk is still sufficiently massive to undergo gravitational 
instability. Under more typical conditions, disk-born objects will likely grow well above the deuterium 
burning planetary mass limit. We conclude that if planets are formed by gravitational instability, 
they must be the low mass tail of the distribution of disk-born companions. To validate this theory, 
on-going direct imaging surveys must find a greater abundance of brown dwarf and M-star compan- 
ions to A-stars. Their absence would suggest planet formation by a different mechanism such as core 
accretion, which is consistent with the debris disks detected in these systems. 
Subject headings: planet formation, accretion disks, binaries 



1. INTRODUCTION 

Motivated by the recent discovery of massive planets 
on wide orbits, we explore the requirements for mak- 
ing gas giants at large separations from their host star 
via gravitational instability, hereafter, GI. In particu- 
lar, we consider the formation mechanism for the system 
HR 8799 which contains three ~10A/jup objects orbit- 
ing at distances between ~ 30 and 70 AU (Marois et al. 
[2008 ) . The standard core accretion model tor planet tor- 
mation, already strained in the outer solar system, has 
difficulty explaining the presence of these objects. While 
GI i s an unlikely f ormation mechanism for close in plan- 
ets ( |Rafikov|2005| , for more widely separated planets, or 
sub- stellar companions, the viability Gl-driven fragmen- 
tation deserves further investigation. 

In the inner regions of a protoplanetary disk, gas can- 
not cool quickly enough to allow a gravitationally unsta- 
ble disk to fragment into planets ( Rafikov|2005 Matzner 
fc Levin|2005l ). For this reason, core accretion— m whicti 



stars have yet to turn up similar companions (Nielsen 



solid planetesimals collide and grow into a massive core 
which then accretes a gaseous envelope — has emerged as 
the preferred mechanism for forming planets at stellar 
separations < 10 AU. Planets at wider separations have 
only recently been disc overed by direct im aging around 
the A-stars HR 8 799 ( |Marois et "al] |2008D , Fomalhaut 
( Kalas et al. 2008), and possibly Beta Pic ( [Lagrange et al. 

Searches at large radii surrounding solar-type 



& Close 2009). Standard core accretion models cannot 
lorm these planets, though further investigation is war- 
ranted. 

In favor of this possibility, all three systems show some 
evidence of processes related to core accretion: all have 
infrared excess due to massive debris disks at large radii. 
This is at least partially a selection effect as these systems 
were targeted due to the disks' presence. Nevertheless, 
these debris disks are composed of reprocessed grains 
from collisions of planetesimals. The disks' long life- 
times prohibit a primordial origin for small grains — they 
are removed quic kly by radiation pressu re and Poynting- 
Roberston drag (Aumann et al. 1984) and so must be 
regenerated from collisions between larger bodies that 
formed through the coagulation of solids at early times. 
Therefore, planetesimal formation, a necessary ingre- 
dient in co re accret ion models, has taken place (e.g., 
Youdin fc S hu 2002| [Chiang fc Youdin[[2009| . 



In addition, other A-stars host planets at 1-2 AU 
([Johnson et al. 2007) which, although they have a dis- 



tinct semi-major axis distribution from planets orbit- 
ing G and M stars, are likely formed by core accretion. 
If future surveys demonstrate that this distribution ex- 
tends smoothly to wide separation planets, then simplic- 
ity would argue against a distinct formation mechanism 
for the wide giants. 
Yet the standard core accretion model faces a serious 



2 



problem at large distances. The obser ved lifetimes of gas 
disks are short, at most a few M yr ( Hillenbrand et al. 
1992 Jayawardhana et al. 2006). Iii contrast, typical 



The planets around HR 8799 pro be a previously un- 



core accretion times inc rease with radius and exceed 10 
Myr beyond 20 AU (e.g.|Levison fc Stewart'2001 ; Goldrc- 



pFh et al.,2004j ) . Whether or not this theoretical ditticuity 
can be overcome will require careful modeling of the in- 
teractions between planetesimals and the young gas disk. 

Could wide orbit plan ets have formed at smaller radii , 
and migrated outwards? Dodson- Robinson et al. (2009), 
have investigated the possibility ot torming the HR 8799 
system via scattering in the absence of dynamically im- 
portant gas, but find that putting three massive plan- 
ets into such close ly spaced yet wide orbits is unlikely. 



circumstances outward migration in resonances might be 
feasible. Alternatively, the core of a giant planet could 
be scattered outward by a planet, or migrate outward 
before accreting its gas envelope, either by i nteractions 
with the gas dis k (Type HI migration; e.g. Masset & 
Papaloizou|2003 ) or with planetesimals embedded in the 



gas (Capobianco, Duncan, & Levison, in prep). Neither 
mechanism has yet been shown to move a core to such 
large distances, though this possibility has not been ruled 
out. 

Given these difficulties, it is natural to search for other 



formation mechanisms, and GI (Boss 1997) stands out 
as a promising alternative. If any planets form by GI, 
the recently discovered directly imaged planets ar e the 

,,2009 iNero 



most likely candid ates ()Rafikov|j2009[ |Boley„200! 
fc Bjorkman|2009[ ). In t nis paper, we examine thii 



is possi- 

bility in more detail, considering the expected mass scale 
of fragments and the effect of global disk evolution on the 
formation process. 

The inferred masses for the HR 8799 plan ets are close 
to the deuter ium burning limit of IBMjup (Chabrier & 
Baraffe|2000 )■ For simplicity, we take this as a the divid- 



ing line between planets and brown dwarfs and we refer 
to the HR 8799 objects as planets throughout. However, 
there is no reason for a given formation mechanism to 
function only above or below this threshold, and in fact, 
we will argue that if the HR 8799 planets formed by 
GI, their histories are more akin to those of higher-mass 
brown dwarfs than to lower-mass planets. 

To constrain GI as a mechanism for wide giant planet 
formation, we set the stage by describing the HR 8799 
system in f|2] We review the standard requirements for 
fragmentation in ^and discuss the initial mass scale of 
fragments in f|4] In S|5] we show that under typical disk 
conditions, fragments will continue to accrete to higher 
masses. We then discuss important global constraints 
on planet formation provided by star formation models, 
disk evolution timescales, and migration mechanisms in 
^ and f|7] We compare predictions of our analysis with 
the known wide substellar companions and exoplanets 
in fpj suggesting that future observations will provide 
a dennitive answer to the formation mechanism for HR 
8799. In the appendices we re-examine the heating and 
cooling properties of disks that are passively and actively 
heated, with special attention to the implications for ir- 
radiated disks, which become increasingly relevant for 
more massive stars. 

2. THE HR 8799 SYSTEM 



explored region of para meter space (Marois et al 



2008 



Lafreniere et al. 2009) because they are more distant 
from their host star. The three companions to HR 
8799 are observed at separations from their host star 
of 24, 38, and 68 AU. Their masses, estimated using 
the observed luminosities of the planets in conjunction 
with cooling models, have nominal values of 10, 10, 
and 7 Mjup, respectively. A range in total mass of 
19-37 Mjup is derived from uncert ainties in the age of 



the host star (Marois et al. 2008). Interpretation of 



Crida et al. ( 2009 ) have suggested that under favorable W 



the cooling models generates substantial additional sys 
tematic uncertainty — recent measurements suggest that 
these mod els may overpredict t he masses of brown dwarfs 
25% ( |Dupuy et al.|2009| ). [Fabrycky fc Murray-Clay 
(120081 have demonstrated tnat tor planetary masses in 
the stated range, orbital stability over the age of the sys- 
tem requires that the planets occupy at least one mean- 
motion resonance, and that for doubly-resonant orbital 
configurations total masses of up to at least 54Mjup can 
be stable. 

HR 8799 has been called a "scaled-up solar system" 
in terms of the stellar flux incident on its giant planets 
( Marois et al.|[2008 Lafreniere et al.|[2009 ). However for 
understanding the lormation ot this system it is more 
useful to consider dynamical times, and disk mass re- 
quirements. Because the dynamical time at fixed radius 

1/2 

scales only as Af* , the dynamical times are larger at the 
locations of the HR 8799 planets that at the solar system 
giants. Since the total mass in planets greatly exceeds 
the ^ 1.5 Afjup in the solar system, we can infer that 
(as with some other extrasolar systems) the primordial 
disk around HR 8799 was more massive than the solar 
nebula and/or there was greater efficiency of planet for- 
mation, especially in the retention of gas. Compared to 
solar system giants, longer dynamical times make core 
accretion more difficult and larger disk masses make GI 
more plausible. 

3. IDEAL CONDITIONS FOR GI-DRIVEN FRAGMENT 
FORMATION 

We first determine where, and under what local disk 
condition s, fragmentation by GI is possible. Following 

and I 



possible. I'qiiqwing 
( |2005|) and [Rafikovl 
th surface density ^ 



,Gammie| (j2001|), |Matzner fc Levinj 

, we argue that tor a disk with surface density U 
and temperature T to fragment, it must satisfy two crite- 
ria. First, it must have enough self-gravity to counteract 
the stabilizing forces of gas pressure and rotational shear, 
as quantified by Toomre's Q: 



Q 



< Qo - 1 



(1) 



(Safronov 1960 Toomre 1964), where Cg — y^kT/ji is 
the isothermal sound speed ot the gas with mean particle 
weight fi = 2.3mH appropriate for a molecular gas, G is 
the gravitational constant, k is the Boltzmann constant, 
and n is the orbital frequency. Equation ([T]) specifies the 
onset of axisymmetric instabil ities in linear theory that 
can g ive rise to bound clumps ( Goldreich fc Lynden-Bell 
19651 p] 



^ In a realistic disk model, clumps likely form within spiral arms 
formed via non-axisymmetric, non-linear instabilities, although the 
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The second criterion that must be satisfied for frag- 
mentation to proceed is the so-called cooHng time cri- 
terion. The heat generated by the release of gravita- 
tional binding energy during the contraction of the frag- 
ment must be radiated away on the orbital timescale so 
that increased g as pressure does not stall further collapse 
( |Gammie|[200l| ). This implies 



tr. 



)Ol 



37Sc^ /(r) 
32(7 - 1) crT4 



(2) 



Here C is a constant of order unity, 7 is the adiabatic 
index of the gas, and a is the Stcfan-Boltzmann constant. 
We take /(r) = 1/t + t (jRafikov 2005) for disk vertical 
optical depth r = kE/2 and gas opacity k. Numerical 
models of collapse in barotropic disks measure the critical 
value C through the inclusion of a loss term u/tcooi in the 
equation for the internal energy, u. Estimates of C r ange 



from ~ 3 — 12, depending on 7 (Rice et al. 2005 ), th e 



numerical implementation of cooling (Clarke et al.|2007 ), 
and the vertical stratification in the disk. We assume 7 = 
7/5, appropriate for molecular hydrogen, and we adopt 
C — i here. Although ( was measured in disks whose 
temperature is controlled by viscous heating, we show in 
Appendix K that the same expression (modulo slightly 
different coefficients) should apply when irradiation sets 
the disk temperature, as will be the case in disks prone 
to fragmentation (see Appendix [Bj) . 

A disk satisfying Toomre's criterion for instability 
(equation [l]) but not the cooling time criterion (equa- 
tion [2]) experiences Gl-driven angular momentum trans- 
port which regulates the surface density of the disk so 
that Q ~ Qo ~ 1 and Q does not reach substantially 
smaller values (c.f. Appendix [A| . We can therefore use 
Toomre's criterion to define a relationship between S and 
T at fragmentation, as a function of period: 



nGQo 



(3) 



where for convenience we define ee {k/ fiy/^{TTGQo)-^ ■ 
We shall hereafter set Qo — I- 

Given equation (|3]), we can rewrite equation (2]) to gen- 
erate a single criterion for fragmentation whicn depends 
on temperature and location: 



nt, 



cool 



c 



ifjt) ^f{r) < 1, 



(4) 



where ft = (3/32)7(7— l)~^fc(^cr()~^. Somewhat coun- 
terintuitively, the critical cooling constraint requires that 
a disk be sufficiently hot to fragment. The value of /(r) 
depends on both f2 and T. Evaluating this criterion relies 
on the disk opacity, which we return to in j]4.2[ 



4. MINIMUM FRAGMENT MASSES AND SEPARATIONS 

4.1. Initial masses of Gl-born fragments 

We take the initial mass of a fragment to be the mass 
enclosed within the radius of the most unstable wave- 
length, Aq = 2ttH in a Q = 1 disk, or: 



frag 



S{2ttH)^ 



(5) 



critical value of Q at which fragmentation occurs should remain 
similar. 



(Levin 2007 ), where iJ = Cs/^ is the disk scaleheight. 



Cossins et al. ( |2009 ) have shown that even when the 
Gl is non-axisymmetric, the most unstable axisymmetric 
wavelength Ag is one of the dominant growing modes, 
suggesting that this is a reasonable estimate. While 
more numerical follow up will be necessary to pin down 
the true distribution of fragments born through GI, at 
present simulations show that this estimate may we ll be 
a lower l i mit, but is the correct order of m agnitude (Bo- 



ley|[2009l iStamatell os fc Whitworth||20 09'). 



smg equation (|3|), we rewrite the fragment mass ex- 
plicitly as a function of temperature and location: 



M, 



frag 



47r 



H 



M* = fr. 



y3/2 



(6) 



where /,„ = {2'k)'^ fqk / fi. Equation ^ demonstrates that 
at a given disk location, fragment masses depend only on 
temperature, with lower temperatures generating smaller 
fragments, subject to the minimum temperature required 
fo r fragmentatio n by equation (H . 



Rafikov (2005) pointed out tnat there exists an abso- 
lute minimum fragment mass at any disk location, when 
the disk satisfies the equalities in equation ([T]) and equa- 
tion ([2]) and has r = 1. The minimum temperature re- 
quired for fragmentation scales as T cx (r-f 1/r)^/^ (equa- 
tionj4]), so the critical temperature and fragment mass are 
minimized at r = 1, the optical depth for which cooling 
is most efficient. The corresponding minimum mass as a 
function of location is: 

Mf,,„i„ = /,„(/, /O'^'f^'/' (7) 

-3/10 / M \ 1/10 



1 



5Mjup (: 



.100 AUy \1.5Mq 
which occurs for disk temperatures: 



^f''"'" = ^^(lOO^) 



-6/5 



1.5Mf; 



2/5 



(9) 



Equation ([7| corresponds to Q = 1, riicooi = 3, and 
T — 1. This minimum mass is only achieved for tem- 
peratures consistent with Tf i^in- Once an opacity law 
is specified which relates T and r, the problem becomes 
overconstrained — these three criteria can only be satis- 
fied at a single disk location, and equations 7]|9 are valid 
at only one radius in the disk. We now proceed to eval- 
uate the critical disk temperatures and fragment masses 
for realistic opacity laws, demonstrating that planetary- 
mass fragments can only form at large separations from 
their host star. 



4.2. Opacity 

At low temperatures, when T < 155 K, ice grains ar e 
the dominant source of opacity (Pollack et alj 1985) 

thi 



Above ~155 K, ices begin to sublimate. In the cold 
regime applicable to the outer regions of protoplanetary 
disks, we consider two realistic opacity laws, one which is 
characteristic of grains in the interstellar medium (ISM), 
and one which is characteristic of grains that have grown 
to larger sizes due to processes within the disk. We as- 
sume a Rosseland mean opacity which scales as 



(10) 



4 



where the both the exponent, /3, and the constant 
are not well constrained in protoplanetary disks. They 
depend on the number-size distribution and composition 
of the dust grains, and on the dust-to-gas ratio. For ISM 
like grains, the Rosseland mean opacity may be approx- 
imated by an opacity law with /3 — 2, or 



(11) 



Pollack et al. 1996 Bell & Linl 1994 Semenov et al. 



is obs ervationally confirmed in the ISM (Beckwith et al. 
2000 and references therein). For our fiducial model we 
use K2 ~ 5 X 10~^c m^/g for T in Ke l vin, a fit to the stan- 
dard dust model by Semenov et al. (2003). Throughout, 
we quote k per gram of gas tor a clust-to-gas ratio of 
10-2. 

As grains grow larger, they eventually exceed the wave- 
length of the incident radiation, and so the opacity is de- 
termined by the geometric optics limit. In this case the 
Rosseland mean opacity is independent of temperature, 
so the exponent (3 — 0. In this limit: 



(12) 



For a fixed dust mass in grains of size s, kq oc s~^. For 
our fiducial model we choose kq sa 0.24 cm^/g which is 
valid for T > 20 K and typical grain sizes of order 300/j,to 
(see Figure 6 of [Pollack et^[l985] ) . 

Observations of emission from optically thin protoplan- 
etary disks show evidence of grain growth at millimeter 
wavelengths. Specifically the measured opacities Kj, oc v"' 
with a ~ 0.5 — 1.5 in the Rayleigh- Jeans tail imply par- 
ticle growth toward the millimete r wavelengths of the 
observations ( D'Alessio et al.|[2001 ). Although most ob- 
served disks have had more time for grain growth to pro- 
ceed, Class sources also show evidence of grain growth 
( D'Alessio et al.|2001 ). Alternatively, these objects could 
be optically thick at rnillimeter wavelengths, which would 
mimic the effects of grain growth. 

Using these opacity laws we see that cooling proceeds 
with a different functional form in the optically thick 
and optically thin limits. In the optically thick regime, 
/(r) T — {npT^Yj) /2, which when combined with 
equation (|4| indicates that to fragment, the disk must 
have temperature 



1/(2-/9) 



1/(2-/3) 



(13) 



for /3 7^ 2. In the special case /3 = 2, the fragmentation 
constraint is temperature-independent, as discussed in 

ine optically thin regime, /(r) = l/r — 2/{KpT^T,). 



Int 



The cooling time is independent of E, so we can rewrite 
equation (|4| as: 



T > {2ftf'^^+l''^ 



^y/(3+/3) 

up) 



(14) 



Fig. [T] shows the dependence of the cooling time on 
disk temperature for each opacity law at two different 
radii. Since fragmentation can only occur when iltcooi < 
C, the minimum temperatures at which fragmentation is 
allowed are specified by the intersection of the cooling 
curves with the fiicooi boundary. 
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OOap as long as the ice grains dommating the opacity 



are smaller than a few tens of microns. This opacity law G 
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Fig. 1. — The disk cooling time as a function of tempertature 
for different opacity laws at radii of 40 AU (dashed) and 100 AU 
(solid). The cooling time is calculated assuming that Q = 1. The 
temperature independent (large grain) opacity law is shown in red, 
while the ISM opacity law: re (X is shown in blue. The line 
thickness indicates the optical depth regime. When lines drop be- 
low the critical cooling time (grey), disk fragmentation can occur. 
The bend in the ISM opacity curve indicates that in the optically 
thick regime, the cooling time becomes constant as a function of 
temperature. 

4.2.1. Small grain opacity law 



For an optically thin disk with (3 = 2, equation ( 14 1 



implies that in order to fragment, the disk must have 
temperatures in excess of: 



r>Tthin =9K 



-3/10 



^2 

5 X lO-^cmV 



100 AU. 

-1/5 



(15) 



M, 
1.5M, 







1/10 



Colder disks, even when optically thin, cannot cool 
quickly enough to fragment. 

For an optically thick disk, /? = 2 turns out to be 
a special case: the temperature dependence drops out 
of equation ( 13 ), giving instead a critical radius beyond 



which fragmentation can occur, independent of the disk 
temperature: 



> 70 AU 



1.5Mf; 



1/3 



K2 



5 X 10-4cm2/g 



-2/9 



(16) 



Matzner fc Levin ( 2005 1 first pointed out the existence 
of a minimum critical radius for fragmentation. In Fig. [2] 
we illustrate how the two fragmentation criteria create 
a radius rather than temperature cutoff. At radii larger 
than the critical radius defined above, any S — T combi- 
nation which gives Q < 1 will fragment (so long as the 
opacity law remains valid). At smaller radii, no combina- 
tion of S and T which gives Q <\ will fragment because 
the disk cannot simultaneously satisfy the cooling time 
criterion. 

4.2.2. Large grain opacity law 
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Fig. 2. — Fragmentation can only occur in the region of param- 
eter space indicated by the overlapping hashed regions for ISM 
opacities at radii of 100 AU. The upper, shaded region (red) shows 
where Toomre's parameter Q < 1. The lower, shaded region (blue) 
indicates where tcool < 3r2~^. At radii less than 70 AU, fragmenta- 
tion is prohibited because the two regions no longer overlap. That 
the boundaries of these regions are parallel lines reflects the k qc 
form of the ice-grain-dominated opacity at low temperatures. 



As shown in Fig. [T] for the large grain opacity, t > 1 
for all relevant temperatures. In this case equation ( 13 1 
requires: 



T > 6bK 



AL r 



1.5M, 



43 AUy V0.24cni2/g 
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The corresponding minimum mass for these tempera- 
tures is: 



P V 1.5M0 

3/4 



5/8 



(43Au) (1^) 



0.24cm2/g 



We have scaled equations (17) and (18) to an effective 
critical radius for this opacity law. Although fragmen- 
tation can occur inside 43 AU at sufficiently high tem- 
peratures, the fragments exceed 13Mjup, making it irrel- 
evant for planet formation. Smaller values of kq move 
this boundary inward, allowing for fragmentation into 
lower mass objects at smaller radii, although the scal- 
ing with opacity is shallow. Grain growth to larger sizes 
could plausibly reduce kq. If, for example, grains domi- 
nating the disk opacity have grown up to 1mm without 
altering the dust-to-gas ratio, equation ( 12 ) implies that 
in the geometric optics limit kq = 0.072cm2/g. In this 
case , the minimum radius is p ushed inward to 26 AU (see 
also [Nero fc Bjorkman||2009 1 . 

Thus tar we have determined the minimum masses al- 
lowed as a function of radius with the temperature as a 
free parameter. We now calculate actual disk tempera- 
tures, which at large radii are typically higher than the 
minima. In this case we must evaluate fragment masses 
using equation ([6]). 



4.3. Initial fragment masses with astrophysical disk 
temperatures 

To consider the case favorable to GI planet formation, 
we consider the lowest plausible disk temperatures in or- 
der to minimize the fragment masses. We estimate the 
disk temperature using the passive flared disk models of 
Chiang & Goldreich (1997). This model likely underesti- 
mates the temperatures m actively accreting systems be- 
cause it ignores significant "back heating" from the infall 
envelope and surroundi ng cloud (Chick & Casscn 1997 
Matzner & Levin 2005). Although viscous heating wil. 
also contribute to the temperature, we ignore its modest 
contribution to obtain the lowest reasonable tempera- 
tures and fragment masses. Disk irradiation dominates 
over viscous heating in this regime (see Appendix B). 

We consider the inner region where the disk is opti- 
cally thick to blackbody radiation. In this regime, the 
temperature of a flared disk in radiative and hydrostatic 
equilibrium is: 



T 

J. u 



r 



1/2 



n OC ^2/7^-3/7 



(19) 



where ap measures the grazing angle at which starlight 
hits the disk; ai? is dependent on the degree o f disk flar 
ing measured at the height of the photosphere ([Chiang 
Goldreich|19 97') . Grain settling may reduce the hejight 
the photosphere (set here to 4 times the scaleheight) and 
thus ap- For the standard radiative equilibrium model, 
the disk flaring scales approximately as H/r (x r^/'^. We 
shall flnd when we calculate the disk temperature that 
a gravitationally unstable disk remains optically thick, 
justifying the use of this formula. 

To estimate the the s tellar luminosity we use the stel - 
lar evolution models of Krumholz & Thompson (2007), 
which include both nuclear burning and accretion energy. 
The accretion luminosity depends on both the current 
accretion rate and the accretion history ( in so far as 
it effects the stellar radius), so we obtain the lowest lu- 
minosity estimates by allowing the star to accrete at a 
constant, low accretion rate. We use the stellar luminos- 
ity after accreting to 90% of its current mass (or 1.35M0 
assuming roughly 10% is still in the disk). We choose an 
accretion rate of lO~^M0/yr as a lower bound because 
a star accretingat a lower accretion rate throughout its 
history has a formation timescale that is too long. Accre- 
tion rates an order of magnitude larger give comparable 
luminosity (when the star has only reached 1.35Mf)) t o 
the present day luminosity of 5Lq ( Marois et al.| |2008). 
Lowering the accretion rate below this value does not 
significantly lower the stellar luminosity because the ac- 
cretion energy contribution is small. 

The luminosity calculated for an accretion rate of 
lO^^M0/yr translates to temperatures of: 



T« 40 



70 AU 



-3/7 



(20) 



which we shall use as our fiducial temperature profile. 
In the outer regions of the disk where fragmentation is 
allowed, the disk temperatures are of order 30 — 50K. 
These temperatures exceed the minimum threshold for 
fragmentation, and so the mass of fragments will be set 
by equation (|6|. 
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Fig. 3. — Depiction of the current configuration of HR 8799 and 
formation constraints for realistic disk temperatures. We show 
the lowest expected irradiated disk temperatures (blue) and cor- 
responding fragment masses (grey), as a function of radius. The 
lower bound on both re gime s (burgundy) is set by the irradia- 



tion model described in [4.3 with M = 10~'^ Mq/jv. The upper 
boundary is set by the current luminosity of HR 8799, ~ 5Lq . The 
green dashed-dotted line shows the mass with disk temperatures 
of 10 K, a lower limit provided by the cloud temperature. The 
vertical line shows the critical fragmentation radius for the ISM 
opacity law; fragmentation at smaller radii requires grain growth. 
Fragment masses are shown for radii at which the irradiation tem- 
peratures are high enough to satisfy the cooling time constraint of 
equation \17\ . At smaller radii, fragmentation is possible at higher 
disk temperatures, but the resulting fragments have correspond- 
ingly higher masses, and planet formation is not possible. 

Other analytic and numerical models of stellar irradi- 
ation predict temper atures in agreement with or higher 
than our estimate. ( [Rafikov fc D e Colle| [2006) jOffner] 
et al.|2 0()9). Similarly, models ot disks in Ophiuchus have 
similar temperatures for 1 Myr old stars of lower mass 
(and thus luminosity), implying that ou r model temper- 
ature s are low, though not unrealistic (Andrews et al. 
2009| 



In Fig. [3] we illustrate the constraints on fragment 
masses from this irradiation model, calculated using 
equation ( ^ We show the fiducial disk temperatures of 
equation ( 5D I , along with temperatures consistent with 
luminosities up to the present-day luminosity. For our 
fiducial opacity laws, the expected fragment masses are 
only marginally consistent with GI planet formation — 
fragments form near the upper mass limit of ISAfjup. 
Lower opacities produced by grain growth might al- 
low fragmentation into smaller objects at closer radii. 
Whether grain-growth has proceeded to this extent in 
such young disks is unclear. 

5. GROWTH OF FRAGMENTS AFTER FORMATION 

For realistic disk temperatures, it is conceivable that 
fragments may be born at several Mjup. We now con- 
sider their subsequent growth, which may increase their 
expected mass by an order of magnitude or more. 

The final mass of a planet depends sensitively on nu- 
merous disk properties (effective viscosity, column den- 
sity, scaleheight) and the mass of the embedded object. 
In order to constrain the mass to which a fragment will 
grow, we can compare it to two relevant mass scales: the 
disk isolation mass and the gap opening mass. 




Fig. 4. — Contours of the ratio of planetary isolation mass to 
stellar mass as a function of Toomre's Q and the disk aspect ratio 
H/r, illustrating that the isolation mass is always large in unsta- 
ble disks. For disks with higher Q's consistent with core accretion 
models, the isolation mass remains small. The shaded region indi- 
cates where the isolation mass exceeds the stellar mass. 

Halting the growth of planetary mass objects is a rele- 
vant problem independent of the formation mechanism. 
However the GI hypothesis requires that the disk is (or 
was recently) sufficiently massive to have Q ~ 1, imply- 
ing that the disk is actively accreting. The core accretion 
scenario does not face this restriction. 

5.1. Isolation Mass 

We estimate an upper mass limit for fragments by as- 
suming that they accrete all of the matter within several 
Hill radii: 

Miso - ^TrfH^Rnr. (21) 

Here fn ^ 3.5 is a numerical constant representing from 
how many Hill radii, Rh = "''(A'^iso/SM*)^/^, the planet 
can accrete (Lissauer||1987 ). 



It is instructive to compare the ratio of the isolation 
mass to the stellar mass: 



Mi, 



H 



3/2 



(22) 



Wc find that large isolation masses are always expected 
in gravitationally unstable disks. Fig. [4] illustrates the 
scaling of equation ( 22 ) with Q and the disk aspect ratio, 
H/r. For our fiducial temperature profile, H/r « 0.09 
at 70 AU. For low values of Q and comparable H/r, the 
isolation mass exceeds 10% of the stellar mass. For the 



ideal disk values cited above (equation 161, the isolation 



Miso « 400Afjup 



70 AU 



6/5 



l.bM, 



o 



1/10 



(23) 



This mass is nearly two orders of magnitude larger than 
the minimum mass. Growth beyond the isolation mass is 
possible either through mergers or introduction of fresh 
material to accrete by planet migration or disk spreading. 

Objects which grow to isolation mass cannot be plan- 
ets, and so we turn to mechanisms that truncate frag- 
ment growth below the isolation mass. 

5.2. Gap opening mass 

Massive objects open gaps in their disks when gravita- 
tional torques are sufficiently strong to clear out nearby 
gas before visco us torques can replenish the region with 
new material. (Lin & Papaloizou 1986 Bryden et al. 



7 



19991. The gap width is set by the balance between the 

(24) 



two torques: 



A 

r 



1/3 



where A i s the gap width, fg ^ 0.2 3 is a geometric factor 
derived in Lin fc Papaloizo u (|1993 ), q is the planet to star 
mass ratio, and a is the SEakura & Sunyaev ( 1973 ) ef- 
fective viscosity. This can be used to derive the standard 
minimum gap opening mass by requiring that A > H: 



/3^ 



(25) 



4^10 (or) Un^ (toau) 



5/4 



1.5 M, 







-5/4 



Gap opening requires A > i?H and Rn > H . The lat- 
ter requirement is automatically satisfied for fragments 
formed by GI. 

While the ef fects of GI are often paramete rized through 
an a viscosity, Balbus k. Papaloizou ( 1999 ) have pointed 
out that a, a purely local quantity, may not adequately 
descri be GI driven transport , which is inherently non- 



local. Lodato & Rice (2005) have shown that for suf- 
ficiently thin disks, the approximation is reasonable: in 
order to form objects of planetary mass, the disk must be 
relatively thin and at least marginally within this limit. 
However, even in this thin-disk limit, it is not clear that 
GI driven torques will exactly mimic viscous ones at gap- 
opening scales 



Equation ( 25 ) implies that the gap opening mass is 
less than or equal to the fragment mass for effective vis- 
cosities consistent with G I. We use a ^0.1, as this is 
consistent with active GI (lGam mie|2001 Lodato & Rice 



2005||Krumholz et al. 2007). If the local disk viscosity is 
lower, fragments will always form above the gap-opening 
mass. 

5.2.1. Gap-opening starvation mass 

Gap-opening slows accretion onto the planet, but does 
not starve it of material completely. Accretion rates 
through gaps remain uncertain for standard core ac- 
cretion models, and numerical models are not available 
for accretion onto the distended objects formed through 
GI fragmentation. Nevertheless simulations of accretion 



through gaps in low viscosity disks (Lubow et al. 19991 
demonstrate that accretion is slower through larger gaps, 
and this qualitative conclusion likely remains valid as 
long as gaps form. 

Analogous to the isolation mass, we consider a "gap 
starvation mass" that is related to the ratio of the gap 
width to planet Hill radius. Rewriting equation (24) we 
find the ratio of gap width to Hill radius is: 



A 



■fgl 



1/3 



(26) 



Note that A > Rh recovers the canonical gap opening 
estimate appropriate for Jupiter: q > 40h'/{r^n), modulo 



order unity coefficients (cf. Crida et aL]|2006 ). 

If we make the simplifying assumption that gap accre- 
tion terminates when the gap width reaches a fixed num- 
ber of Hill radii, fs, we can calculate a gap starvation 
mass. We expect that for a gap to truncate accretion, 
fs ^ fn = 3.5, the width in Hill radii used to calculate 
the isolation mass ( |5.1[ ). Terminating accretion is an un- 
solved problem for Jupiter in our own solar system, and 
so to provide a further c onstraint on fs , we re fer to the 
numerical simu lations of Lissauer et al. (2009) (See also 



D'Angelo et al. (2003) for a detailed explanation of the 



numerical work)." Iheir runs 21 and 21 J exhibit asymp 
totic mass growth after ^^2.5 Myr for a constant-mass, 
low viscosity (a = 4 x 10~^) disk under conditions appro- 
priate to the formation of Jupiter. Using equation (26), 
we solve for the width of the gap generating this fall-off 
in accretion rate and find fs = A/Rh ^ 5. The need for 
an extremely large and well-cleared gap reflects the inte- 
grated effects of low-level accretion through the gap and 
onto the planet over the disk lifetime of a few Myr. Even 
a slow trickle of material onto the planet can contribute 
to significant growth. 

Using /s = 5, the gap-opening starvation mass for HR 
8799 scaled to both the simulated solar-system viscosity 
and to the expected GI viscosity is: 



'-8M 



J up 



40 K / V70 AU/ 



4 X 10- 



/ A 



\5Rh 



(27) 



^2000M,„,(^)(^)' 

(—) ( — )■ 
\40 K J V70 AUy 



In order to limit growth to planetary masses, the effec- 
tive viscosity must be two orders of magnitude below 
that expected in GI unstable disks, roughly a ^ 10^^. 
More restrictively this requires that other local trans- 
port mechanisms such as the MRI be weaker than cur- 
rently predicted by simulations — they pro duce a ^ 10 
at least in disks with a net magnetic flux ([Fleming et a^ 



2000 Fromang et aL]|2007 Johansen et al.||2009 



5. Fig. I 

tn radius 



illustrates the scaling of gap starvation mass wit 
for several values of a. It appears that active disks face 
severe obstacles in producing planetary mass objects un- 
less the disk disappears promptly after their formation. 
Although we expect fragmentation to make the disk more 
stable by lowering the local column density, there is little 
reason to expect a recently massive disk to be so quies- 
cent. 

5.2.2. Planet starvation through gap overlap 

Although it is unlikely that the disk will fragment into 
a sufficiently large number of plan etary mass objects to 
completely deplete the disk of mass ( Stamatellos & Whit- 
[worth ,2009 ) , the formation of multiple fragments simul- 
taneously may limit accretion through competition for 
disk material by opening overlapping gaps. The cur- 
rent separation between the planets is such that gaps 
larger than roughly three Hill radii overlap, so depend- 
ing on their migration history, this could limit growth 
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1000 




1994 1. Fragmentation at earlier times leads to the forma- 



Radius (AU) 

Fig. 5. — The gap starvation mass as a function of disk radius. 
We show curves for several values for a, and indicate the planetary 
mass regime, and the region in which disk fragmentation is likely. 
We use fc: = 5. sc aled to simulation 21J of Jupiter formation in 
[Lissauer et al.| | [2009| ) (labeled L09 in the figure). The radial scaling 
is derived assuming H/r oc r'^^'^ . For the low viscosity case, wo 
normalize the scale height to Jupiter at 5.2 AU in a 115 K disk for 
comparison with L09. For the higher viscosities, we normalize the 
dis k scale height to the lowest expected temperatures (equation 
\20\ . For comparison we show the HR 8799 planets as black circles. 

(see also Q. From equation (27 1 we see that if gaps 
are forced to be smaller by a factor of two due to com- 
petition with another planet, the expected masses are 
decreased by a factor of 8. This effect would imply that 
fragments in multiple systems should be lower in mass. 
Note that when simulated disks fragment into multiple 
objects simultaneously, they generally have orbital con- 
figurations like hierarchical multiples rather than plan- 
etary system s (Stamatellos & Whitworth 2009| Kratter 



et al 
limit 



2010) 



Whether the same conditions required to 
ragment growth — reduced disk mass and/or vis- 
cosity after fragmentation — can allow the retention of a 
planetary system of fragments remains to be simulated. 

5.3. Disk dispersal as a means to limit fragment growth 

Mechanisms for gas dispersal such as photoevapora- 
tion may be necessary to stunt planetary gro wth, even 
for models of Jupiter in our own solar system ( Lissauer 
et al.| 2009 ) . Dissipation timescales for A- sta r gas disks 
are thought to be short, less than 2-3 Myr (Carpenter 



et al. 2006), which could halt growth before the gap 
opening starvat ion mass is achieved. Radiative transfer 



models such as Gorti & Hollenbach (2009) have calcu- 
lated that photoevaporation by the central star will be- 
come important at radii of ^ 100 AU aro und one Myr for 
an A star (see also Ercolano et al."2009). This timescale 
coincides with the expected fragmentation epoch, and 
may aid in shutting off accretion both onto the disk, and 
onto the planets. 

6. GI PLANET FORMATION IN THE CONTEXT OF STAR 
FORMATION 

We now consider how the disk can reach the fragmenta- 
tion conditions described in f|3] in the context of a model 
for star formation. Due to the effects of infall onto the 
disk, we find that planet formation via GI can only occur 
when the fragmentation epoch is concurrent with the end 
of the main accretion phase, as the protostar t ransitions 
from a Class I to a Class II object ( Andre fc MontmerIe| 



tion of more massive companions, while fragmentation at 
later times is unlikely because disks are too low in mass 
( [Andrews et~aL]|2009[ ). 



6.1. Ongoing accretion and the formation of binaries 
and multiples 

Because the cooling time constraint is easily satisfied 
for expected disk temperatures, disks are likely driven to 
fragmentation by lowering Q. Even when Q is above the 
threshold for fragmentation, torques generated by self- 
gravity (e.g. spiral arms) can drive accretion. When 
the infall rate onto the disk is low, a self-gravitating 
disk regulates its surface density and hence Q so that 
the torques are just large enough to transport the sup- 
plied mass down to the star, thereby avoiding fragmen- 
tation. However, GI cannot process matter arbitrarily 
quickly because the torques saturate. Thus the disk will 
be driven toward fragmentation if the infall rate becomes 
too high. This critical accretion rate is a function of disk 
temperature: 

M„it « (28) 



( Gammie 2001 Matzner fc Levin 2005 ). Numeri cal simu- 
lations show that GI saturates at ctsat 0.3— 1 (Gammie 



200T|JKrumholz et al.|2007l|Lodato fc Rice|2005l [Kratter 



et al. 



2010) 



If the infall rate onto the disk exceeds Merit 
the disk can no longer regulate the surface density to 
keep Q just above unity, and fragmentation will occur. 
Because the conversion of accretion energy to thermal 
energy at large radii is inefficient, disks cannot restabi- 
lize through heating to a, rbitrarily high accretion rates 
([Kratter & Matznerj|2006[). 



It fragmentation occurs due to rapid accretion, as de- 
scribed above, it is difficult t o limit subsequent fragment 
growth. As demonstrated in §5.2. 1[ gap opening does not 
limit accretion efficiently when the effective viscosity, in 
this case, ctsat, is high. 

A more general barrier to making small fragments dur- 
ing rapid infall is the large reservoir of mater ial passing 
by the frag ment as star formation proceeds ( Bonnellfc] 
Bate 1994k. The specific angular momentuni, j, of ac- 
creting disk material typically increases with time (mod- 
ulo small random fluctuations in a turbulent core) , land- 
ing at a circularization radius, rdrc = j'^/GM^,, which 
is larger than the fragment's orbit. Since the fragment's 
Hill radius is roughly 10% of the disk radius, newly ac- 
creted material undergoes many orbits in the fragment's 
sphere of influence as it tries to accrete onto the central 
star, and some fraction will accrete onto the fragment 
itself. This process is less efficient if global GI modes 
drive fragments to smaller radii, because their Hill radii 
shrink. However, growth to stellar (or sub-stellar) masses 
may still occur as long as migration timescales are not 
faster than accretion timescales. The latter scenario im- 
plies that the disk cannot fragment at early times and 
still reproduce a single-star system like HR 8799. 

The trend toward continued growth and even mass 
equalization following disk fragmentation is obs erved in 
numerous simulations wit h ongoing accretion (|Bonnell 



fc Bate[[19 94^ 'Bate"2000! 'Matsumoto fc Hanawal |iSOOi}[ 
Walch et al.. .2009^ Kratter et al.,,2010, .Krumholz eFliL 
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2009 1 . Simulations of planet formation by GI with ong o- 
mg accretion also illustrate this behavior ( Boley|[2009 ) . 

Consequently, if the HR 8799 planets formed by GI, 
they must have fragmented at the tail end of accretion 
from the protostellar core onto the disk. Most likely, this 
requires that the protostellar core have properties such 
that its infall rate reaches the critical value in equation 
( 28 ) just as the core is drained of material. If this coinci- 
dence in timing does not occur, fragmentation produces a 
substellar, rather than a planetary, companion. Whether 
fragmentation at this epoch can produce non-heirarchical 
orbits like HR 8799 remains to be explored. 

6.2. Reaching Instability in the absence of infall: FU 
Orionis outbursts 

Driving the disk unstable with external accretion cor- 
responds to excessive growth of fragments. It is there- 
fore tempting to consider mechanisms to lower Q through 
disk cooling, while holding the column density fixed. Be- 
cause the disks are dominated by irradiation, changes 
in the viscous dissipation due to accretion are unlikely 
to affect a significant temperature change, and so lower- 
ing Q requires order of magnitude changes in the stellar 
luminosity due to the weak scaling of T (x p^^"^ ■ FU 
Orionis type outbursts (Hartma nn fc Kenyon||1996 1 can 
cause rapid changes in luminosity, lb result in planetary 
mass fragments, the luminosity drop following an FU Ori 
outburst would need to reach at least the minimum lumi- 



nosity used in equation ( 20 1 on timescales shorter than 



an outer disk dynamical time 

The accretion of GI formed embryos onto the protostar 
is a p roposed source of the outbursts (Vorobyov & Basu 
20061. If this occurs, perhaps a final generation of grav- 
itational fragme nts, the so-called "last of the Mohicans" 
( Gonzalez 1997 ) , would remain as detectable compan- 
ions. Although the lack of infall in this scenario might 
ease the gap opening constraints, fragments face all of 
the other difficulties discussed above in remaining low in 
mass. 

7. MIGRATION IN A MULTI-PLANET SYSTEM 

A final consideration for making wide orbit planets 
through GI is their subsequent migration history. If 
formed via GI, the HR 8799 planets had to migrate in- 
ward to their current locations. There is an independent 
reason to believe that migration did in fact take place 
in th i s syst em. As discussed by [Fabry cky fc Murray-] 
Clayl ( |2008D , the long term stability of HR 8799 requires 
a resonant orbital configuration, most plausibly a 4:2:1 
mean motion resonance, which likely resulted from con- 
vergent migration in the protoplanetary gas disk. We 
demonstrate below that while this history is plausible, it 
requires special disk conditions. 

Inward Type II migration (appropriate for gap opening 
planets) is expected at formation, because the fragment- 
ing region is within the part of the disk accreting onto 
the star. Type II migration of a single planet occurs on 
the disk viscous timescale: 



13/14 



M* 



1.5Mr. 



1/2 



(29) 

(^)"i») 



where we have used v — acj./fl. This timescale is short 
enough to allow substantial migration during the lifetime 
of the gas disk. 

If fragments migrated inward independently, they 
could never become captured in resonance, as the inner- 
most planet would migrate farther and farther from its 
neighbor. However, as shown in 35] planetary gaps may 
overlap if they are within several Hill radii of each other, 
comparable to the current separations between the HR 
8799 planets. This overlap alters the torques felt b y, and 
thus the migration of, the planets. As shown by |Kley| 
([2000 ), if multiple planetary gaps interact, convergent 
migration is possible. Gap interaction allows an outer 
planet to shield an inner planet from the material, and 
thus torques, of the outer disk, slowing or halting its in- 
ward migration and allowing the outer planet to ca tch 
up. This mechanism is invoked by |Lee fc Peale (2002) to 
generate convergent migration and resonance capture in 
the planets orbiting GJ 876. 

Under the assumption that gap overlap allows conver- 
gent migration and resonance capture, we now ask: What 
is the overall direction of the subsequent migration? We 
note that if the planets migrated a substantial distance 
after resonance capture, eccen tricity damping by the gas 
disk was likely necessary (c.f., Lee fc Peale|2002| . We do 
not consider eccentricity damping further here. Once two 
planets are caught in mean-motion resonance, the torque 
on an individual planet from the gas disk can cause both 
planets to migrate, with angular momentum transfer me- 
diat ed by the resonan ce. ^Masset fc Sncllgrove^(2001 ) (see 
also Crida et al. 2009 ) have argued that the torque imbal- 
ance on~a^aIrorgap-opening resonant planets can even 
reverse the direction of migration, although this relies on 
a significant difference between planet masses. 

Nevertheless, understanding the planets' overall migra- 
tion requires understanding how overlapping gaps alter 
the torque balance on the group of planets. Guided by 
our interest in clean gaps that limit the growth of planets 
(^HHP' consider the following simplified problem. 

We imagine that the three planets have cleared, and 
are embedded in, a single large gap which is sufficiently 
clean that any disk gas passing through is dynamically 
unimportant. Because the system is locked in a double 
mean-motion resonance, we assume that an imbalance in 
the torques acting on the two edges of the gap can cause 
all three planets to migrate. We can then ask: is there a 
sufficient flux of angular momentum from the outer disk 
to cause the planets to migrate inward with the viscous 
accretion of the disk? 

When in the 4:2:1 resonance, the total angular mo- 
mentum of the planets is roughly 2.4Afpf2prp, where we 
have assumed that the planets are roughly equal in mass, 
Mp, while Tp and Hp are the separation and Keplerian 
angular velocity of the outermost planet. The angular 
momentum flux from the outer disk is large enough to 
move the planets on a viscous time Tp/i/ when: 



(31) 



where M — SttYiV is the mass flux through radius r^. 

Using the above inequality, we can calculate a critical 
disk surface density at the current location of the outer 
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planet such that the disk can push the planets inward: 



of different masses because disk f ragmentation becom es 



p 



(32) 



For Mp ~ lOAfjup and our fiducial disk temperatures 



(equation 20 1 , this constraint is always satisfied when 
Q = 1. I'he disk is unable to cause inward migration 
when S < 4g/cm^ at 70 AU, which is equivalent to Q ^ 
20. 

If the planets do share a clean common gap, a large 
fraction of the disk would be effectively cleared of gas 
while a massive outer disk is still present. A simi- 
lar mechanism has been invoked to explain transitional 
disks, which contain holes at radii of a few tens of AU 



and smaller (Calvet et al. 2002). Transport of disk gas 
through a less well-cleared gap could substantially alter 
this picture. 

In summary, it is possible to envision a scenario in 
which the HR 8799 planets migrate inward to their cur- 
rent locations in such a way that their orbits converge, 
allowing resonance capture. This scenario is consistent 
with other constraints on GI planet formation: shortly 
after formation, the disk must have low accretion rates 
and decline in mass in order to (a) limit the growth of 
fragments, (b) allow for large, overlapping gaps. 

More stringent constraints will require future work on 
migration in gravitationally unstable disks, particularly 
in the presence of multiple planets massive enough to 
clear large, overlapping gaps. 

8. CURRENT OBSERVATIONAL CONSTRAINTS 

While there is a regime of parameter space in which 
planet formation is possible, typical conditions produce 
more massive (> ISMjup ) companions. If GI fragmenta- 
tion ever forms planets, then fragmentation should typ- 
ically form more massive objects. Consequently, these 
planets would constitute the low mass tail of a distribu- 
tion of disk-born companions. If the mass distribution 
is continuous there should be more sub-stellar compan- 
ions than planets at comparable distances of 50-150 AU. 
Observing this population is a strong constraint on the 
formation mechanism, but current data are insufficient 
to draw conclusions. 



Zuckerman & Song (2009) have compiled the known 
sub-stellar companions m this range of radii to date. This 
range of separations falls be yond the well-established 
inner brown dwarf "desert" (McCarthy & Zuckerman 



2004 1 , and has not been well probed due to observational 
difficulties at these low mass ratios. Note that the overall 
dearth of brown dwarf c ompanions to solar mass stars 
is not a selection effect (Metchev & Hillenbrand 2009 



Zuckerman fc Song||2009[ ) 

We illustrate the ob servational constraints by p lotting 
the companions from Zuckerman & Song (2009) along 
with the known exoplanets compiled by the Exoplanet 
Encyclopedi£(^ as a function of mass ratio and projected 
separation in Fig.lG^, and as a function of minimum frag- 
ment masses and fragmentation radii in Fig. [6b. We com- 
pare these with the HR 8799 planets, Fomalhaut b, and 
the solar system giants. We distinguish between stars 

^ October 2009, [http://exoplanet.e u compiled by Jean Schnei- 



der 



more likely for higher mass stars ( Kratter et alT][2008 1 . 

At present, neither the population of substellar com- 
panions nor the population of exoplanets is continuous 
out to HR 8799. Many selection biases are reflected in 
Fig. |6j and these gaps in particular may be due to se- 
lection effects; resolution and sensitivity make it diffi- 
cult to detect both wide orbit planets, and close-in low 
mass brown dwarfs. We note that while there are actu- 
ally fewer planets at distances less than 1 AU, the cutoff 
above 5 AU is unphysical. There is not yet any indication 
of an outer cut off radius in the exoplanets: if they con- 
tinued out to larger separation, the distribution would 
easily encompass the HR 8799 and Fomalhaut systems. 

Data from ongoing surveys like that which found HR 
8799 are necessary to verify the true companion distri- 
bution as a function of mass and radius. If these planets 
are formed through GI, we would expect observations 
to fill in the gap between HR 8799 and higher mass ra- 
tio objects to show a continuous distribution. If these 
planets are formed via core accretion, than observations 
may fill in the plot on the opposite side of HR 8799, oc- 
cupying a region of parameter space for which neither 
core accretion nor GI is currently a successful formation 
mechanism. 

9. SUMMARY 

We have demonstrated that while GTdriven fragmen- 
tation is possible at wide distances from A stars, frag- 
ment masses typically exceed the deuterium burning 
"planet" limit. In contrast, the formation of sub-stellar 
and stellar companions is more likely because moderate 
disk temperatures and active accretion onto and through 
the disk drive disk-born objects to higher masses. 

If the HR 8799 planets did form by GI, the following 
criteria had to be met: 

1. Fragments should form beyond 40 — 70 AU: inside 
of this location the disk will not fragment into plan- 
etary mass objects even if Q ^ 1. Grain growth is 
required for fragmentation at the lower end of this 
range. 

2. Temperatures must be colder than those of typical 
disks to limit the initial fragment masses. 

3. The disk must be driven unstable at a special time: 
infall onto the disk must be low, but the disk must 
remain massive (e.g. the end of the Class I phase). 
The disk must only become unstable to fragmenta- 
tion at this point because earlier episodes of insta- 
bility should lead to sub-stellar or stellar compan- 
ion formation. 

4. The subsequent growth of fragments must be lim- 
ited through efficient gap clearing necessitating low 
disk viscosity or early gap overlap. Disk dispersal 
via photoevaporation may also be necessary. 

5. The three fragments must form at the same epoch 
separated by several Hill radii, implying that the 
entire outer disk becomes unstable simultaneously. 

6. Migration must be convergent. This likely requires 
the gaps of the planets to overlap so as to starve 
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Fig. 6. — (Left) Known substellar companions (stars) and planets (plusses) as a function of mass ratio and projected separation. The 
three objects in the HR 8799 system are shown by pink circles, and a pink triangle denotes the upper-limit mass ratio for Fomalhaut b based 
on the dynamical mass estimate of'Chiang et al. (20091). Grey squares indicate the gap regions. Ongoing suryeys are necessary to determine 
whether there is a continuous distribution between Jupiter/Saturn (blue diamonds) and HR 8799, or HR 8799 and brown dwarf companions. 
Planets around very low mass primaries with A/, = 0.02-0. IMq are marked by purple squares. These systems haye mass ratios more akin 
to the substellar companions than to the remainder of the population of planets. Primary masses range from M« = O.O2M0-1M0 (black) 
and M« = IMq— 2.9Mq (red) for substellar companions and from M* = 0.1— 0.4Mq (purple) and Af« = 0.4— 4.5AfQ (black) for pla nets . 
(Right) The same objects plotted as function of the minimum fragment mass, Affrag,min a-nd critical radius, rcrit- We use equation (|16| 
to calculate r^rit- For Affrag.mim we apply equation (|6| at radius r^rit under the simplified assumption that the disk temperature is set by 
the stellar luminosity: L/Lq = {Mj, / Mq)^-^ for M, > 0.43AfQ and L oc M^-^ for lower-mass stars. The temperatures used to calculate 
fragment masses are not allowed to dip below 20K. Masses below Affrag,min ^'I'e unlikely to result from GI. 



the inner most planet of disk material, thereby pre- 
venting runaway inward migration. 

If these conditions are met, then the planets in HR 
8799 could comprise the low-mass tail of the disk-born 
binary distribution, the runts of the litter. In this case 
one would expect to find a larger number of brown dwarfs 
or even M stars in the same regime of parameter space - 
surrounding A-stars at distances of 50 — 150 AU. 

Ongoing direct imaging surveys of A and F stars will 
provide a strong constraint on the formation mechanism 
for this system: if HR 8799 is the most massive of a 
new distribution of widely separated planets, our anal- 
ysis suggests that formation by GI is unlikely. On the 
contrary, the discovery of a population of brown dwarf 
and M-star companions to A-stars would corroborate for- 



mation via disk fragmentation. 
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APPENDIX 



A. COOLING AND FRAGMENTATION IN IRRADIATED DISKS 



Rafikov ("2009) has suggested that the cooHng time might be ahered in an irradiated disk. Here we consider the 
cooling time for thermal perturbations to a disk, and show that a simple formula [equation ( A7 1] gives the cooling time 



for arbitrary levels of irradiation, at least in radiative optically thick disks. 

Consider ambient radiation striking an optically thick disk with a normal flux Fq — aT^ = (3/8)-Firr (where the 
numerical factor in the last definition is purely for later convenience). Note that To depends both on the irradiation 
field (from the host star and/or light emitted and reflected from a surrounding envelope) and also on the disk's 
surface geometry, e.g. flaring angle. By incorporating these variables into Tq we attempt a general calculation. At the 
photosphere, where the optical depth to the disk's self-emission r = Tphot ~ li energy balance gives: 



aT^f ^ aT^ + F/2 . 



(Al) 



Here, F is the luminous flux from any internal sources of energy dissipation, e.g. viscous accretion, shocks, or the 
gravitational binding energy released by a collapsing fragment. The factor of two reflects that half of the radiation 
is emitted from each surface of the disk. Since op tical light is absorbed abo ve the IR photosphere, about half the 
irradiation free streams out before it heats the disk ( Chiang fc Goldreich|1997 |. We can absorb this reduction into the 
definition of T^. 
We assume heat is transferred by radiative diffusion as: 



4 dT^ 

ri7 



F 
2"' 



(A2) 



since convection is suppressed by irradiation and may be a negligible correction in any event (Rafikov 2007). Integration 
from the midplane at r = Ttot and T = to the photosphere gives: 



(4/3)a(r4 - T^g) = (F/2)(rtot - Tphot) ■ 



(A3) 
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We now drop the subscripts from rtot and the midplane Tm, which we wih soon take as the char acteristic temperature 
(ignoring order unity corrections from height averaging). Furthermore we apply equation (Al) and take the Tphot ^ 
Ttot — ^ T hmit (meaning that we don't need to know the precise location of the photosphere] to express 



8a 
37 



(T* - 



1 /8ctT4 



F. 



(A4) 



which shows that the midplane temperature is controlled by the larger of Ft and i^irr- 
The cooling timescale to radiate away thermal fluctuations (generated e.g. by GI) is 

_ T.6U 
^cool - -j^ 



(A5) 



where a temperature perturbation ST has an excess heat SU « cpST, and cp = {k/ /[■^ — 1) is the specific heat for 
a mean molecular weight /i and adiabatic index 7. Strongly compressive motions, which arc not at constant pressure, 
will introduce order unity corrections that we ig nore. 
The excess luminous flux, using equation (A4 1, is 



(5F = 



8ct 5T 
37 7^ 

32crr3(5r 

37 



[(4-/3)r4 + /3r4] 



(l-/?/4) if r»To,/3^4 
1 if T ~ To 



(A6) 



where t = kS/2 cx . The point is that the escaping flux varies by only an order unity factor between the strongly 
(T ~ To) and weakly (T 3> To) irradiated regimes. Typical grain opacities, < /? < 2, ensure the correction is order 
unity (and also ensure that we ca n ign ore the catastrophic heating that would occur if /3 > 4). 
Combining equations ( A5 ) and ( A6 1 with the definition of heat capacity we find that the cooling time is simply 



''COOl 



37Sc2 



32(7 - l)aT^ 



(l-/3/4)-i 
1 



if T»To,/?^4 
if T ~ To 



(A7) 



where the isothermal sound speed Cs = \JkT j [i. 

From this derivation, we see that the cooling time obeys the simple form of equation ( |A7[ ) for all levels of irradiation 
— which only introduces an order unity /? correction. The cooling time depends on pTor weakly irradiated disks 
because changes to the opacity alter the amount of flux that escapes from the midplane. In highly irradiated disks, 
opacity changes have little effect because the small difference between the midplane and surface temperatures drives 
a weak flux. 

While increasing the irradiative flux incident on a disk decreases the cooling time by raising T, it will not trigger 
fragmentation in a Q ^ 1 disk, since it increases Q. We will not explore optically thin or convective disks at this time. 

We assume that iltcooi < C ~ 3 is the fragmentation criterion inde pendent of irrad iation. When the cooling time 
is longer, the disk is presumed to enter a state of gravito-turbulence (Gammie 2001). We take this term to mean a 
quasi-steady state of gravitationally driven turbulence, on scales < n wherein viscous dissipation of GI turbulence 
regulates Q ~ 1. Thus the cooling time can be tra nslated to a the critical value of a at which gravito-turbulent 
accretion disks will fragment. While jRice et al. (2005) flnd that an a-threshold is more robust than one for ^cooi when 
the adiabatic index varies, they did not include irradiation, which we contend would reveal that cooling is ultimately 
the more physical criterion, but the issue is best settled by simulation. The emitted flux with an a-viscosity and Q ^ \ 
gives 

,3o2 in (A8) 



F K. (9/4)i^Sf7^ w 9/(47r)ac^f77G 



which combined with equation ( A4 ) and equation (|A7| without the /3 correction gives 



f7t 



cool 



7 



7- ly 4(F + Ti,,/T) V7-1 



7 



1 

9a 



1 



F- 

f7 



(A9) 



When irradiation is weak enough that F\„ < Tr, we recover the standard a ^ 1 criterion for fragmentation (ignoring 
the accumulated order unity coefficients). However for stronger irradiation with i^irr ^ Ft^ fragmentation occurs for 
a > Ft IF\^^^ a lower threshold. 

Gravito-turbulent models require modification when Fi„ > F ^ i.e. an even lower level of irradiation than needed to 
affect the a fragmentation threshold. In this case the disk shows some similarities to isothermal disks, and should have 
lower amplitude density perturbations, because there i s in sufficient viscous dissipation to support order unity thermal 
perturbations (see the related discussion in|Rafikov 2009]) . We note that simulations of iso thermal disks do develop 
GI, exhibit Gl-driven transport and fragment ( jKrumhoiz et al. 2007 Kratter et al. 20101, but they do not appear 
particularly turbulent. 
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B. TEMPERATURE DUE TO VISCOUS HEATING 



We now show that visc ous heating is relatively unimportant in the outer reaches of irradiated A-star disks (see also 
section 3 of Rafikov|2009 1 . For an optically thick disk with an ISM opacity law, balancing viscous heating and emitted 
radiation gives 



— aT^ « Am/22 (b^) 
3r Stt 



The solution for the midplane temperature is 

2/3 / • \ 2/3 



1287r2 aGQo I \ IQ-SMo/yr / ^° V70 AU/ 

lower than the irradiation temperatures shown in Fig. ISlby more than a factor of four. Note that the surface density 
falloff 

is also quite steep for a constant Qo ^ 1, viscous disk with the ISM opacity law. 
If the disk is optically thin, then the balance between heating and cooling gives: 

AraT^^-^AW^ (B4) 
on 



i9 K 



lO-6Af0/yr 



This temperature profile is shallow, but still colder than the irradiation temperature at large radii where the disk 
becomes optically thin. 



